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This  report  covers  research  which  was  carried  out  from  October  1,  1978 


through  September  30,  1979  under  Grant  No.  AFOSR  77-3158.  The  results 
are  in  four  main  areas:  (1)  Frequency  domain  tests  for  optimality  in  problems 
of  periodic  control,  (2)  Infinite -horizon  problems  as  a  natural  extension  of 
traditional  problems  in  perioidic  control,  (3)  Applications  of  periodic  control 
to  aircraft  cruise,  (4)  Minimal  realization  of  nonlinear  (functional  2-power) 
input/output  maps.  The  attached  bibliography  includes:  items  which  have 
been  reported  in  previous  Interim  Scientific  Reports  and  are  required  for 
reference  here  [1,  2,  5]  ,  new  items  [3,  6]  ,  and  reports  currently  in  prepara¬ 
tion  [4,  7]  .  The  significance  and  interrelationship  of  these  items  and  other 
new  results  is  reviewed  in  the  following  paragraphs.  Elmer  G.  Gilbert  was 
the  principal  investigator.  Daniel  J.  Lyons  and  Dennis  S.  Bernstein,  doc¬ 
toral  students  at  The  University  of  Michigan,  also  worked  under  the  Grant 
and  made  important  contributions. 

The  v  test  [a]  is  a  second  order  test  for  optimality  in  periodic  control 
problems.  It  involves  sinusoidal  perturbations  in  the  neighborhood  of  a  steady- 
state  optimum  and  provides  a  powerful  algebraic  test  which  may  determine 
whether  or  not  a  system  is  proper  (time  dependent  periodic  control  gives  bet¬ 
ter  performance  than  steady-state  control).  Reports  [1,2]  extended  the  work 
in  [a]  in  three  directions:  (1)  applied  the  ir  test  to  a  more  general  periodic 
control  problem,  (2)  gave  conditions  which  assured  the  validity  of  the  n  test 
([a]  contained  crucial  omissions),  (3)  explored  in  detail  relationships  between 


second  order  conditions  for  optimality  in^l^e j^^e^dyj-jSftate  and, dynamic  problems. 
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During  the  last  year  the  material  has  been  extended  further  in  two  ways:  the 
gap  between  necessary  conditions  and  sufficient  conditions  for  optimality  has 
been  closed  (this  was  based  on  [b]  ,  whose  applicability  was  discovered  last 
fall),  the  proofs  (contained  in  [2]  )  have  been  simplified.  The  extended  re¬ 
sults  are  now  contained  in  a  single  report  [3]  which  will  appear  as  a  regular 
paper  in  the  December,  1979  issue  of  the  IEEE  Transactions  on  Automatic 
Control. 

The  preceding  work  assumes  that  the  control  input  is  unconstrained 
(u(t)  €  R111).  A  result  which  applies  in  the  presence  of  constraints 
(u(t)  €  U  c  R3 n)  has  been  obtained  this  year.  It  relies  on  a  basic  lemma  in 
a  recent  paper  by  Warga  [c]  and  provides  a  test  for  proper.  It  can  be  sum¬ 
marized  somewhat  imprecisely  by  some  additions  to  the  notation  of  [3]  . 
Introduce 


jt  =  *  (0)  +  (CG(0)  +  D)'  Hyy(CG(0)  +  D), 

$  +  Tt  («  ))  (£  -  V  (W» 

jt(u)  = 

(jr  -  it  («))  (ff  +  tt(w)) 

Then  the  periodic  control  problem  is  proper  if  there  exists  a  pair  of  vectors 
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This  "generalized  ar  test"  goes  beyond  the  rather  special  case  treated  in 
[d]  .  Several  examples  have  been  examined  which  demonstrate  that  it  is  a 
useful  test. 

The  new  test  requires  a  normality  condition  which  is  more  complex 
than  the  condition  given  in  [3]  .  Considerable  effort  was  given  to  finding  a 
more  easily  verified  condition.  It  proved  only  partially  successful.  The 
present  theory  requires  U  to  be  convex,  but  it  is  clear  that  this  condition 
may  be  relaxed  in  a  variety  of  ways.  The  study  of  these  and  other  related 
questions  is  now  nearing  completion. 

In  the  usual  periodic  control  problem  the  cost  and  constraints  depend 
on  a  vector  y  (see  [3]  for  notation)  whose  components  are  given  by 

y.  =  ~  J  f.(x(t),  u(t))dt,  i=l»  •  •  • » m, 

1  0  1 

where  x(t),  u{t)  are  (periodic)  solutions  of 
x(t)  =  f(x(t),  u(t)),  x(0)  =  x(t). 

When  the  period  t  is  large  there  is  some  question  as  to  whether  or  not  the 
formula  for  y^  is  valid  and  the  motion  should  be  periodic.  An  alternative 
formulation  is 

oo  A  -6  .t 

Yi  ~  S  f  Wt),u(t))e  1  dt,  6  >  0, 

0  1 

”6it 

where  e  represents  the  exponential  discounting  of  future  costs  to  the  present. 
Some  economists  have  considered  problem  formulations  in  this  direction,  but 
apparently  without  knowledge  of  results  in  periodic  control  (see,  e. g.,  [e]  ). 
The  theory  of  the  above  "infinite  horizon"  problem  has  many  similarities  with 
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the  theory  of  periodic  control.  In  particular,  there  is  a  steady-state  optimi¬ 
zation  problem  contained  within  the  dynamic  optimization  problem:  x(t),  u(t) 
constant  gives 

f(x, u)  =  0  and  y.,  =  (6.)  1  ^(x,  u). 

Thus  it  is  possible  to  compare  steady-state  performance  with  dynamic  per¬ 
formance.  Preliminary  investigations  have  produced  some  interesting  results. 

For  example,  if  x(t)  and  u(t)  are  periodic  with  period  t  (certainly  a 
special  choice  for  the  infinite  horizon  problem),  it  can  be  shown  that 

-6  .T  T  -6  .t 

y.  =  (1  -  e  *  )’  J  e  *  f  (x(t),  u(t))dt. 

For  t6^  <<  1,  this  gives 

-  t  . 

y.  =  T  J  £  (x(t),  u(t))dt. 

1  0 

Thus,  the  usual  periodic  control  problem  can  be  viewed  as  an  infinite  horizon 
problem  in  which  the  motion  is  periodic,  the  period  is  small  (6^t  «  1),  and 

~  - 1  A 

f  j  =  (6.)  f..  Alternatively,  the  tools  of  periodic  control  theory  may  be  ap¬ 

plied  to  the  study  of  the  infinite  horizon  problem.  This  has  been  done  with  a 
variation  of  the  n  test  and  necessary  conditions  for  the  dynamic  optimality 
of  steady -state  solutions  have  been  obtained.  This  technique  is  also  useful 
in  examining  steady-state  solutions  different  than  the  one  suggested  above, 
e. g. ,  those  which  appear  when  x(0)  is  fixed  (see  [e]  ).  Many  open  questions 
remain  and  will  be  the  subject  of  research  in  the  coming  year. 

The  computational  work  on  periodic  cruise  of  aircraft  was  continued. 

The  aircraft  model  used  is  described  in  [f]  .  Previous  computations 
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concerned  the  optimization  of  specific  range  and  assumed  that  a  (engine  spe¬ 
cific  fuel  consumption)  was  constant.  This  year  a  more  realistic  model  for 
<r  was  considered.  Its  dependence  on  speed  and  altitude  was  a  good  repre¬ 
sentation  of  a  typical  small  jet  engine.  The  resulting  optimal  periodic  mo¬ 
tions  were  close  to  those  obtained  for  <r  =  constant.  With  fairly  strong  con¬ 
straints  on  the  maximum  altitude  (say  20,  000  ft.  or  less)  the  improvements 
in  performance  were  notably  better  than  for  a  =  constant.  The  poor  steady- 
state  cruise  performance  of  the  realistic  engine  at  the  constraint  altitude 
was  the  principal  factor  in  this  difference. 

In  addition,  the  problem  of  maximizing  endurance  in  cruise  was  inves¬ 
tigated.  The  aircraft  model  was  the  one  described  in  [f]  with  <r  =  constant. 
Analysis  of  the  energy-state  approximation  gives  a  good  indication  of  the 
principal  factors  which  make  periodic  cruise  more  efficient  than  steady-state 
cruise.  It  shows  that  optimum  periodic  cruise  tends  toward  low  altitudes 
where  maximum  engine  thrust  is  greater  and  flight  speed  (drag)  is  lower. 
Thus,  it  is  necessary  to  impose  an  explicit  constraint  on  minimum  altitude 
(h  >  h^).  The  figure  shows  the  reduction  in  average  fuel  rate  obtained  by 
periodic  motion:  .TfJ^)  *  is  the  ratio  of  average  fuel  rate  with  periodic 
control  to  the  optimal  fuel  rate  with  steady-state  cruise,  Tm(Tgs)  *  is 

the  ratio  of  maximum  engine  thrust  to  the  thrust  required  for  maximum 

2  -1 

endurance  cruise,  P  =  V#(gho)  where  is  the  steady-state  speed  for 
maximum  endurance.  Optimum  periodic  performance  using  the  full  point- 
mass  model  gives  the  points  designated  by  +  and  x.  Several  computations 
have  been  made  using  the  realistic  engine  model.  They  show  smaller 
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improvements  than  those  indicated  in  the  figure.  However,  if  an  altitude 
constraint  of  the  form  hQ  <  h  <  h^  *8  imposed,  a  much  greater  improvement 
over  the  optimum  steady -state  cruise  is  obtained.  This  suggests  that 
periodic  cruise  can  greatly  extend  endurance  when  flight  close  to  the  earth' s 
surface  is  required. 

The  research  on  periodic  control  of  aircraft  cruise  is  largely  due  to 
D.  J.  Lyons  and  will  appear  in  [4]  ,  which  is  in  final  stages  of  preparation. 
This  report  should  lead  to  several  articles  in  the  open  literature. 

Investigations  in  nonlinear  systems  theory  have  focussed  on  minimal 
order  realizations  for  2  -power  input/output  maps  of  the  form 
oo 

y(t)  =  /  k(T J,  T2)  u(t-T1)u(t-T2)dT1dT2. 

A  solution  to  this  problem  was  presented  in  [5]  for  the  case  where  k  is 
continuous.  Two  extensions  were  obtained  this  year:  the  inclusion  of  im¬ 
pulsive  terms  in  k,  the  treatment  of  the  corresponding  discrete -time  prob¬ 
lem.  The  minimal  realizations  obtained  in  these  extensions  have  a  partic¬ 
ularly  simple  form.  Do  all  other  minimal  realizations  have  a  relation  to 
these  simple  realizations?  This  question  was  answered  affirmatively  (for 
k  continuous)  in  [6]  .  Moreover,  the  ideas  extend  to  impulsive  k  and  the 
discrete-time  case.  The  consequence  is  a  complete  theory  for  the  minimal 
realization  of  2-power  input/ output  maps.  A  paper  which  includes  all  of 
these  results  [7]  is  in  preparation. 
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